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Space weather, including solar storms, can impact Earth by disturbing the geomagnetic 
field. Despite the known dependence of birds and other animals on geomagnetic cues for 
successful seasonal migrations, the potential effects of space weather on organisms that 
use Earth’s magnetic field for navigation have received little study. We tested whether 
space weather geomagnetic disturbances are associated with disruptions to bird migra- 
tion at a macroecological scale. We leveraged long-term radar data to characterize the 
nightly migration dynamics of the nocturnally migrating North American avifauna over 
22 y. We then used concurrent magnetometer data to develop a local magnetic distur- 
bance index associated with each radar station (ABmax), facilitating spatiotemporally 
explicit analyses of the relationship between migration and geomagnetic disturbance. 
After controlling for effects of atmospheric weather and spatiotemporal patterns, we 
found a 9 to 17% decrease in migration intensity in both spring and fall during severe 
space weather events. During fall migration, we also found evidence for decreases in 
effort flying against the wind, which may represent a depression of active navigation 
such that birds drift more with the wind during geomagnetic disturbances. Effort flying 
against the wind in the fall was most reduced under both overcast conditions and high 
geomagnetic disturbance, suggesting that a combination of obscured celestial cues and 
magnetic disturbance may disrupt navigation. Collectively, our results provide evidence 
for community-wide avifaunal responses to geomagnetic disturbances driven by space 
weather during nocturnal migration. 
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Earth’s magnetic field is occasionally but regularly impacted by bursts of energy from the 
Sun (“space weather"), such as coronal mass ejections (1, 2) (Fig. 14). On Earth, the 
impact of larger geomagnetic disturbances caused by space weather includes the auroras 
as well as geomagnetic storms that disrupt satellite communications, navigation systems, 
and power grids (2, 3). Despite the magnitude of interference that large geomagnetic 
disturbances can have on human-built technology, less is known about how geomagnetic 
disturbances affect biological systems, including animals that depend on Earth’s magnetic 
field for migratory orientation and navigation. 

That animals use Earth’s magnetic fields for orientation and navigation has been deci- 
sively shown through decades of observation and experimentation (5, 6). There is strong 
evidence that birds, sea turtles, and other organisms depend on magnetic information at 
multiple spatial scales, keying into small changes in inclination, intensity, and declination 
when making orientation decisions and developing navigational maps (7—10). Experimental 
evidence that birds alter navigation decisions based on minor changes in the local magnetic 
field (7, 11) suggests that large space weather events could potentially disrupt the ability 
of migrating animals to navigate. However, evidence of impacts from naturally occurring 
disturbances on free-roaming animals is scarce. Previous studies in migratory birds have 
suggested potential correlations between geomagnetic disturbances and changes in dis- 
persion patterns of bird headings (12), changes in migratory routes (13) or arrival patterns 
after oceanic crossing (14), shifts in pigeon homing direction (15, 16), increases in bird 
vagrancy (17), and decreased expression of migratory restlessness in captivity (18). Other 
studies have reported relationships between sunspot counts—a correlate of increased solar 
activity and space weather—and both cetacean beach strandings (19) and pigeon homing 
efficiency (20). However, such studies have been relatively localized in geographic, tem- 
poral, or phylogenetic scale, often focusing on specific migration observatories or single 
species (but see ref. (17)). Additionally, some previous results have been inconsistent or 
contradictory, such as those using bird headings and homing pigeons (21, 22), highlighting 
how little is known of how animal navigation is affected by geomagnetic disturbances, 
especially at landscape and population scales. 

Here, we leverage long-term ecological and geomagnetic datasets at a macroecological 
scale to test how nocturnally migrating birds in a continental flyway interact with the 
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dynamism of the magnetosphere. Weather surveillance radars 
capture community-wide dynamics of nocturnal migration, pro- 
viding detailed measurements of migration intensity and direction 
over the United State’ NEXRAD network (23-26). We use a 23-y 
(1995 to 2018) dataset on bird migration collected from 37 radar 
stations in the central migration flyway of the US Great Plains, a 
major migratory corridor (Fig. 1B). We specifically chose this rel- 
atively flat region to minimize variation in migratory orientation 
caused by mountainous topography or oceanic and Great Lakes 
coastlines (27, 28). The community of nocturnally migrating birds 
in this region is primarily composed of a diverse set of perching 
birds (Passeriformes; 73% of species), shorebirds (Charadriiformes; 
12%), and waterfowl (Anseriformes; 9%) (29). We then matched 
data from each radar station with a customized, spatiotemporally 
explicit index of geomagnetic disturbance, ABmax, which repre- 
sents the maximum hourly change in nanoTeslas (nT) from back- 
ground magnetic conditions (Materials and Methods). We derived 
ABmax from time series of geomagnetic measurements collected 
from magnetometer stations near the focal radar stations (Fig. 1B 
and SI Appendix, Figs. S1 and S2 and Table S1) (30). These local 
data, which we accessed from the superMAG database (31), are 
processed to remove daily and seasonal fluctuations, leaving only 
variations caused by space weather effects (Materials and Methods). 
As the dynamics of nocturnal bird migration are greatly affected 
by atmospheric weather conditions (26, 32), we model the impact 
of geomagnetic disturbances on bird migration while controlling 
for the expected effects of atmospheric weather (Materials and 
Methods). 

We hypothesized that geomagnetic disturbances hinder effective 
orientation and navigation during migration. Specifically, we 


predicted decreases in nocturnal migration intensity—the quantity 
of migrating birds—during high geomagnetic disturbance as a 
consequence of increased hesitancy to initiate or continue migra- 
tion with unstable navigation cues. We also predicted decreased 
effort flying against the wind, hypothesizing that birds migrating 
under strong magnetic disturbances will spend less effort actively 
navigating in flight and consequently fly in greater alignment with 
the wind. We further predicted that any effects of geomagnetic 
disturbances on migratory parameters should be magnified when 
other navigation cues are obscured, specifically under overcast 
conditions, which can obstruct nocturnally migrating birds’ abil- 
ities to see celestial cues (33, 34). Finally, a previous study found 
evidence for decreases in the altitude at which migrating birds fly 
in response to a geographic magnetic anomaly, potentially to be 
closer to terrestrial landmarks or to follow changing magnetic field 
parameters (35), so we additionally tested whether similar 
decreases in mean altitude occurred in response to space weather 
geomagnetic disturbances. 


Results and Discussion 


Modeling Framework. We used two complementary methods to 
quantify the influence of geomagnetic disturbances on migration 
response variables while accounting for complex spatial and 
temporal autocorrelation with atmospheric weather predictors. 
First, our “global” nonlinear mixed effects (NLME) models 
include all predictors and produce CI simulated from model error 
parameters (36, 37). These models successfully recovered previously 
documented relationships between migration response variables 
and atmospheric weather and spatiotemporal predictors (23, 32) 
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Fig. 1. Conceptual and geographic layout of our study system. (A) Space weather from the Sun, such as coronal mass ejections, disturbs Earth's magnetic field, 
causing the auroras and potentially decreasing the magnetic field's reliability for migrating birds. Artwork by John Megahan. (B) Distribution of NEXRAD radar 
stations (dark blue circles) and SuperMAG inventory magnetometer stations (purple crosses) used in this study in relation to topography (grayscale) (4). We 
used the three closest and active magnetometer stations surrounding each radar station (S/ Appendix, Fig. S1) to interpolate ABmax, or maximum change in the 
magnetic field from quiet conditions, every hour. Some magnetometer stations had periods of missing data (S/ Appendix, Fig. S2 and Table S1), so we sampled 
magnetometers from a larger geographical area than the radar stations to achieve a robust time series. 
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(SI Appendix, Figs. S3-S8), providing confidence in their predictive 
abilities in this complex dataset. Our second statistical framework 
uses machine learning decision tree (MLT) models to residualize 
our response variables against the weather and spatiotemporal 
variables (38). MLT models have been previously used to produce 
accurate predictions of migration in response to these variables (26), 
and thus, residuals reflect variance not explained by weather and 
spatiotemporal variables. The MLT residualization process excludes 
ABmax, and in models testing an interaction between ABmax 
and cloud cover, it also excludes cloud cover. We then tested for an 
effect of ABmax and cloud cover on the residuals using nonlinear 
regression and covariance matrix-estimated confidence bands (39, 
40). Full modeling details are provided in Materials and Methods. 

We assessed the support for model trends by comparing CI at 
the median ABmax value for each season (noted as ABmax; „Ts 
the same for both seasons) and the highest ABmax value for a 
particular measurement (ranging from 1,690 to 1,870 nT, depend- 
ing on the occurrence of geomagnetic disturbances under varying 
circumstances, noted as e.g., ABmax,¢o9 nr). We interpret non- 
overlapping and narrow CI at these points as showing support for 
a trend, and we verified these CI through visualization of the 
model predictions. NLME comparisons are not centered around 
0 because they result from conditional effect predictions, so 
changes in bird behavior at high ABmax are relative to predictions 
from low ABmax. As such, both numbers may be positive, even 
if the model predicts decreases. 


Migration Intensity and Geomagnetic Disturbance. Our models 
generally recovered support for a relationship between large 
geomagnetic disturbances (ABmax) and decreased migration 


intensity. Migration intensity decreased at large-magnitude 
geomagnetic disturbances by 11.2% in our spring MLT residuals 
model (Fig. 2A and SI Appendix, Fig. S9; 95% CI ABmaxy, „ri 
-0.001, 0.000, ABmax,754 nr: -0.19, -0.03). Consistent with this 
trend, our spring NLME model predicted a 17.4% reduction in 
migration intensity at high ABmax (Fig. 2B), although support 
was slightly weaker in the NLME model as CI at the median 
and highest ABmax overlap partially (Fig. 2B; 95% CI ABmax,, 
at 2-88, 3.106, ABmaxy794 vri 2.68, 2.96). In the fall, our MLT 
residuals model predicted an 8.9% decrease in migration intensity 
at high ABmax (Fig. 2C and SI Appendix, Fig. S10; 95% CI 
ABmax,, „t -0.0003, 0.001, ABmax973 gr -0.14, -0.03). 
These three models predicted decreases in migration intensity 
at geomagnetic disturbances greater than ~500 nT (Fig. 2). By 
contrast, our fall NLME model had low support for a relationship 
between migration intensity and ABmax (Fig. 2D). 

The decreases in migration intensity predicted by our models 
suggest that birds are less likely to initiate or continue migration 
under unstable geomagnetic conditions. Such a relationship sug- 
gests that geomagnetic disturbances may make the magnetic field 
less reliable for birds, affecting their decision to migrate. This idea 
is consistent with experimental simulations of geomagnetic dis- 
turbances that found decreases in migratory restlessness in some 
species (18). However, our models with an interaction between 
ABmax and cloud cover did not consistently recover changes in 
migration intensity at high ABmax and overcast conditions in 
either season (SI Appendix, Figs. S11-S14), suggesting that geo- 
magnetic disturbances do not have a stronger effect on the overall 
number of migrating birds under overcast conditions, contrary to 
our predictions. 
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Fig. 2. Migration intensity predictions with geomagnetic disturbances. (A) Our MLT residual model and (B) NLME model predict decreases of 11% and 17%, 
respectively, in migration intensity at high geomagnetic disturbance (ABmax) in the spring. (C) Our MLT residual model predicts migration intensity decreases 
of 9% at high ABmax in the fall, but the corresponding NLME model (D) does not recover a similarly strong relationship. The spring models and the fall MLT 
residual model begin predicting decreased migration intensity with geomagnetic disturbances of around 500 nT. 
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Effort Flying against the Wind and Geomagnetic Disturbance. To 
test whether birds’ active navigation decreased during geomagnetic 
disturbances, we assessed the impact of geomagnetic disturbances 
on the crosswind component of bird airspeed. Here, the crosswind 
component is the vector component of bird flight, or the flight 
energy, that is perpendicular to the wind vector and represents 
effort flying against the wind. We predicted that birds migrating 
under large geomagnetic disturbances would not navigate as 
effectively and would drift more with the wind, decreasing their 
crosswind component. We found the main effect of ABmax to be 
associated with a decrease of the crosswind component only in the 
fall MLT residuals model (S7 Appendix, Fig. S15). This relationship 
was best supported at high but not extreme ABmax values 
(SI Appendix, Fig. S15; predicted decrease of 1.7%; ABmax;o00 nT: 
-0.030, -0.004). Our models did not find changes in the crosswind 
component in response to ABmax in the fall NLME model 
(SI Appendix, Fig. S16) or in either spring model (SI Appendix, 
Figs. S15 and S16), weakening support for a relationship between 
a main effect of ABmax and the crosswind component. Yet, when 
we tested for an interaction between cloud cover and ABmax, our 
fall NLME model predicted a 25.4% decrease in the crosswind 
component at 100% cloud cover and high ABmax (Fig. 38B; 
95% CI ABmax,; „r: 0.48, 0.54; ABmax,6o4 nr: 0.16, 0.34). 
The MLT residuals model also showed a 7.6% decrease under 
those circumstances, but this relationship has lower support from 
model CI (SI Appendix, Fig. S17; 95% CI ABmax,, „r: -0.007, 
-0.002; ABmaxy¢94 nr: -0.215, 0.054). Evidence for decreases 
in effort flying against the wind with overcast conditions and 
large geomagnetic disturbances during fall is consistent with our 
hypothesis that migrating birds have difficulty navigating during 
geomagnetic disturbances while other navigation cues are obscured 
(33, 34). That is, birds aloft during geomagnetic disturbances 
in the fall may end up drifting with the wind more often. This 
interpretation aligns with research showing that birds drift with 
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the wind upon sudden release at migratory altitudes under cloudy 
conditions, whereas they fly in the expected migratory direction 
against the wind under clearer conditions (41). 

Interestingly, we only found relationships between ABmax and 
the crosswind component in the fall, whereas in the spring none of 
the models showed well-supported trends (Fig. 3A and SI Appendix, 
Fig. S18). Notably, a recent study on the relationship between geo- 
magnetic disturbances and bird vagrancy—individuals appearing 
outside of their expected seasonal range—also found effects only in 
the fall (17). These patterns could be linked if birds migrating under 
conditions faulty for navigation are thrown off course by the wind 
and become vagrants. Our radar measurements of fall migration 
include many inexperienced juveniles migrating for the first time, 
whereas spring migration is composed of birds that have completed 
at least one leg of a migratory round trip. Therefore, the differences 
in spring versus fall patterns shown here may be a consequence of a 
more disruptive effect of magnetic disturbance on the migration of 
inexperienced, young birds (7) which are thought to depend on 
simpler “clock and compass” magnetic orientation strategies rather 
than the “true” magnetic map navigation displayed by experienced 
migrants (5, 42-44). There may be other nonexclusive explanations 
for these seasonal differences, all of which require further research. 
The well-documented faster pace of spring versus fall migration, 
which is thought to be driven by competition for nesting sites and 
phenological demands (45—47), may select for greater effort flying 
against the wind in spring than fall, regardless of magnetic condi- 
tions. Additionally, prior work proposes that birds may allow for 
more drift if they are far from their destination and compensate for 
drift as they approach their goal (48). This idea is in line with the 
seasonal differences that we detect here because in spring long-distance 
Neotropical migrants observed within our study region are likely 
closer to their destinations than they are in fall. Finally, at the mid- 
latitudes that encompass our study area, large geomagnetic storms 
cause local decreases in the magnetic field strength (49), although 
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Fig. 3. Crosswind component of airspeed predictions with geomagnetic disturbances at different cloud cover levels. (A) Our spring NLME model does not 
predict changes in the crosswind component with magnetic disturbance or different cloud cover levels. (B) Our fall NLME model predicts a 25% decrease in 
effort flying against the wind at high ABmax and 100% cloud cover. The model otherwise predicts no changes in the crosswind component during geomagnetic 


disturbances at lower cloud covers. 
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exceptions occur due to auroral and other currents (50). When the 
local magnetic field strength decreases, it could temporarily mimic 
the magnetic field from further south, potentially causing complex 
seasonal differences in how birds respond to geomagnetic 
disturbances. 


Flight Altitude and Geomagnetic Disturbance. We found little 
support for an effect of geomagnetic disturbances on the mean 
altitude of migrants aloft. In both seasons, models that included 
only a direct effect of ABmax on mean altitude did not suggest a 
general relationship between migration altitude and geomagnetic 
disturbances (SI Appendix, Figs. S19 and S20). Models with an 
interaction between cloud cover and ABmax also did not broadly 
suggest interacting effects under overcast conditions (SJ Appendix, 
Figs. S21-S24). These results contrast with decreases in altitude 
observed over a different magnetic stimulus, a geographic 
magnetic anomaly (35), which may have enabled recognition of 
terrestrial landmarks or have been a product of birds following 
shifts in inclination isoclines. The nature of these disturbances 
is also different, as the magnetic anomaly was associated with 
an increase in magnetic field strength (35) while geomagnetic 
storms generally involve decreases in magnetic field strength (49). 
Our results suggest that navigating birds may respond differently 
to geographic magnetic anomalies versus temporal disturbances 
caused by space weather. 


Comparison between Atmospheric Weather and Space Weather. 
Our NLME models allowed us to estimate the effects on response 
variables of each weather and spatiotemporal predictor and therefore 
compare their relative impact with geomagnetic disturbances. 
Consistent with our understanding of migration biology, the biggest 
environmental impacts on migration intensity that we detected 
were associated with atmospheric weather and spatiotemporal 
variables (23, 32). For example, with all other variables held 
equal, warm temperatures in the spring were associated with 
a 45.2% increase in migration intensity compared to the mean 
temperatures (SI Appendix, Fig. S4; 95% CI çc: 2.88, 3.105 99 c: 
3.33, 3.55) and north winds in the fall were associated with an 
around 43.1% increase in migration intensity compared to south 
winds (SI Appendix, Fig. $3; 95% CI 20 mis: 3-38, 3.573 20 mist 2-95, 
3.14). Despite the smaller magnitude of the effects of ABmax that 
we detect on bird migration compared to atmospheric weather, 
log-likelihood tests of NLME models supported the inclusion 
of ABmax and an interaction with cloud cover, with only three 
exceptions, even when CI indicated low support for predicted trends 


(SI Appendix, Table S2). 


Conclusions. Our study provides correlational evidence for a 
relationship between nocturnal bird migration dynamics and 
geomagnetic disturbances. Importantly, our analyses explicitly 
controlled for the effects of atmospheric weather instead of 
examining geomagnetic disturbance in isolation, which has been 
a limitation of previous studies (12, 13, 15-17), (but see ref. 
14). We found broad support that migration intensity decreases 
under high geomagnetic disturbance. However, we found that 
an interaction between cloud cover and geomagnetic disturbance 
did not affect migration intensity and was apparently only 
important for one variable, the crosswind component in the fall. 
These results suggest that the decision to migrate or not during 
geomagnetic disturbances is not generally affected by the ability 
to see celestial cues, but that navigation might be impaired under 
those conditions if birds decide to migrate. 

Our analyses further suggest that correlations between migra- 
tion dynamics and geomagnetic disturbance are mostly evident 
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with the largest geomagnetic disturbances in our dataset, specifi- 
cally those that had variations larger than ~500 nT (e.g., Figs. 2 
and 3). These large disturbances occurred on 81 (3.32%) unique 
nights in the fall and 50 (2.28%) unique nights in the spring in 
our 23-y dataset. Incidences of this magnitude are thought to be 
caused by space weather and solar activity (51, 52). Therefore, 
although the frequency of such events is relatively low, our results 
nevertheless strengthen the link between extraterrestrial impacts 
on the Earths magnetic field and the behavior of nocturnally 
migrating birds. Indeed, 75.5% of the large disturbances included 
in our dataset are directly linked to registered geomagnetic storms 
and all are associated with high planetary magnetic activity 
(SI Appendix, Table S3) (53). 

Our findings highlight how animal decisions are dependent on 
environmental conditions—including those that we as humans 
cannot perceive, such as geomagnetic disturbances—and that these 
behaviors influence population-level patterns of animal movement. 
Our finding that geomagnetic disturbances might impact both the 
number of birds migrating and their directionality strengthens 
previous suggestions of the effects of geomagnetic disturbances on 
bird migration (12, 17, 18). By leveraging the massive amount of 
data collected by radar and magnetometer stations, we were able 
to detect and quantify larger-scale patterns than previously possi- 
ble. Finer-scale patterns, such as species-specific responses or inter- 
actions with life history traits, cannot be studied with our current 
dataset but remain a promising topic of future research, perhaps 
in combination with citizen science projects or bird banding efforts 
(17, 29). Traits that vary across birds that migrate at night through 
our study area (e.g., Passeriformes, Anseriformes) include body 
size, migration distance, and development of migration routes 
through social learning versus genetic inheritance (54), all of which 
could affect susceptibility to geomagnetic disturbances. Our work 
further implies that longer-term impacts to Earth’s magnetic field 
could also affect bird migration, such as the 11-y solar cycles that 
affect the frequency of geomagnetic disturbances (55). Further 
development of migration monitoring technology will allow for 
expansions in our understanding of how space weather and changes 
in the magnetosphere interact with biological systems. 


Materials and Methods 


Nocturnal Migration. To characterize nocturnal migration, we sampled NEXRAD 
scans that are made every 5 to 10 min to extract measures approximately every 
30 min, and we used scans beginning after sunset and ending before sunrise 
(average fall 22.53 scans/night, spring 20.44 scans/night). Data were down- 
loaded from the Amazon Web Services repository (https://s3.amazonaws.com/ 
noaa-nexrad-level2/index.html). We used data from spring and fall migration 
periods from the nights of 1 March to 1 June and 1 August to 15 November (26). 
For each radar sampling event, we built vertical profiles of migration intensity 
(from reflectivity) and migrant speed and direction (from radial velocity) at 100 m 
intervals between 100 m and 3,000 m above ground level (56). We used data out 
to a 100-km radius in the construction of these profiles and elevation scans below 
5.0°. For these profiles of migrant behavior, we removed precipitation contamina- 
tion using the MistNet algorithm (57). MistNet is a convolutional neural network 
that classifies individual sampling volumes as biological or nonbiological. 

We used three response variables that measure the mean behavior of birds 
flying over radar stations. We characterized migration intensity from radar reflec- 
tivity (n, cm?/km?), which measures bioscatter and correlates strongly with the 
number of birds migrating, that is, the intensity of bird migration (58). Mean 
altitude or elevation (km) is the distance above the ground at which most birds 
are detected in the air column. The crosswind component of airspeed (m/s) is the 
vector component of the bird airspeed vector positioned perpendicular to the 
wind vector. It quantifies the extent to which a bird is flying against the wind, with 
smaller values resulting from greater alignment with the wind. To calculate this, 
we multiplied the airspeed vector by the sine of the angle between the airspeed 
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vector and the wind vector. If this angle was larger than 90°, we subtracted it from 
180° before taking the sine. 

To simplify statistical analyses, we summed values of migration intensity across 
the vertical profile for each radar scan and then took the base-10 logarithm. For 
altitude, we used mean altitude weighted by migration intensity across 100-m 
vertical intervals. For the crosswind component, we used mean crosswind com- 
ponent weighted by migration intensity across 100-m vertical intervals for each 
radar scan and took the logarithm after adding 0.001 as a correction for obser- 
vations for which the crosswind component was 0 m/s. Prior to summarizing 
our data, we excluded height bins with missing measurements for migratory 
intensity and bird direction. We excluded 2,202 radar scans in which more than 
16 of 30 vertical intervals were missing data. After summarizing our data, we 
excluded 251 radar scans with 0 migration intensity, which were missing data 
for the crosswind component and altitude, and no RMSE, a measurement of error 
in the direction of migration across height bins. We also excluded two outlier 
observations in which bioscatter was moving at >200 m/s against the wind, 
which is biologically unrealistic and likely erroneous, as the maximum crosswind 
component is otherwise 78 m/s. 


Magnetic Disturbance Index. Our hourly magnetic disturbance index, ABmax, 
measures the maximum change in the magnetic field in nanoTeslas (nT) relative 
to an expected baseline. This index is created from minute-resolution observa- 
tions of the Earth's field as recorded at a set of observatories across North America 
(30). The raw observatory values are processed by the SuperMAG data center to 
remove quiet values (31), including the Earth's intrinsic field contribution and con- 
tributions from Solar Quiet (Sq) current systems in the ionosphere. The remaining 
signal, referred to widely as AB, represents magnetic perturbations from quiet 
conditions due to space weather effects. To create a value indicative of magnetic 
activity over longer periods, we took a rolling hourly maximum of the total vector 
magnitude of the disturbed field, resulting in a magnetometer-specific ABmax. 
Next, we linearly interpolated values to each radar observatory using Delaunay 
triangulation (59), yielding the final, station-specific index. The nearest magneto- 
meter was always within 13° longitude and 12° latitude of the target radar station 
(SI Appendix, Fig. S1). We used all North American magnetometers available on 
SuperMAG between central Canada and Panama (Fig. 18). Not all magnetometers 
were active throughout our entire study period and ABmax could not always be 
interpolated for all radar scans (S/ Appendix, Fig. S2 and Table $1). In the end, 
we removed 6,188 and 3,494 radar scans in the fall and spring respectively with 
missing ABmax values from the final dataset. Our final datasets have 1,725,415 
radar scans in the fall and 1,392,826 radar scans in the spring. 

Our methodology allows for better characterization of local magnetic distur- 
bances than standard global geomagnetic activity indices such as Kp, Dst, or the 
sunspot number (60-63). Ground magnetic perturbations can be highly localized, 
especially at higher latitudes (64) such that regional and global indices often miss 
local perturbations. The 3-hourly interval of the Kp index and low-latitude spec- 
ificity of Dst can also miss mid-latitude signatures of substorms, which are local- 
ized to Earth's nightside with timescales of 3 h or less (65). Our index, ABmax, 
accounts for time and space localization of perturbations as best as possible. The 
occurrence of geomagnetic storms is closely correlated with the ~11-y solar cycle, 
which is defined by the cyclical appearance of sunspots (1). As a predictor, A Bmax 
has a finer temporal resolution and allows for better temporal detrending than the 
sunspot number, which is measured daily and follows a cyclical trend. The sunspot 
number is additionally correlated with general solar activity and has been used in 
studies measuring the impact of radiofrequency noise on magnetoreception (17, 
19). However, solar radiofrequency noise is thought to mostly impact the daylit 
side of the Earth (66) and the frequencies known to disrupt magnetoreception 
(67) might be filtered by Earth's ionosphere (68), which may decrease relevance 
for nocturnal migration. 


Weather Data. We used east-west winds (u wind vector component; m/s), 
north-south winds (v wind vector component; m/s), air temperature (°C), pres- 
sure at mean sea-level (Pa), accumulated precipitation (kg/m’), visibility (m), 
relative humidity (%), and total cloud cover (%) from the North American Regional 
Reanalysis (NARR) (69). Collectively, these variables are known to accurately 
predict bird migration behavior as measured with radar (26). Weather data are 
estimated for 3-hourly intervals in a 32-km grid, so each radar observation is 
associated with the closest NARR timestamp and grid block. Wind measurements 
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(u wind, v wind) and air temperature were taken for each height bin in the radar 
air column, so we used the average across height bins. For the calculation of bird 
airspeed, the two vector components of the wind vector were subtracted from 
the bird groundspeed vector components (also in m/s), yielding bird airspeed 
vector components. 


Modeling Framework. We used two complementary statistical frameworks to 
test our hypotheses. Our models control for the known effects of weather (above), 
temporal variables (specifically, decile of night, ordinal date, and year), and geo- 
graphic variables (latitude and longitude) to measure putative effects of magnetic 
disturbances on bird migration. Our crosswind component and altitude models 
also include migration intensity as a predictor (with a cubic root transformation, 
instead of a log transformation, to avoid negative values), as these measures 
change with larger volumes of birds. Our models also use random effects and 
grouping variables to control for the temporal and spatial autocorrelation in our 
data (see below). All analyses were done in R version 4.2.0 (70). All predictor 
variables were centered and standardized using the “scale” function (70). 


NLME Models. We built our mixed-effect models using the "Imer" function (pack- 
age Ime4) (36) and included nonlinear splines using the "bs" function (package 
"splines”) (70). To build our models, we tested various random effect structures 
(see below) and included fixed effects of all weather, geographic and temporal 
variables. We tested random effects on models with an initial fixed effect structure 
in which all predictors had splines of three degrees of freedom except for time of 
night and day of year (which both had splines with five degrees of freedom) and 
latitude, longitude, and year (which had no splines, because they do not vary as 
much between observations). All models were fit with maximum likelihood to 
facilitate model comparison. 

We tested the influence of the following random effects on model perfor- 
mance based on likely sources of autocorrelation in our data: 1) A "nightly" term 
that grouped each night of sampling from each radar station derived from the 
interaction of ordinal date, year, and radar station. This term included a random 
slope with a spline along time of night, both with and without a random intercept. 
This was our most important source of autocorrelation, as we expect measure- 
ments at each radar station for any given night to depend the most on each other. 
2)Anightly term shared across radar stations from the interaction of ordinal date 
and year that also included a random slope with a spline along time of night, 
both with and without a random intercept. This accounts for larger sources of 
autocorrelation (e.g., weather) that might be shared across multiple radar sta- 
tions on any given night. 3) A “seasonal” term from the interaction of year and 
radar station with a random slope and spline along ordinal date, also both with 
and without a random intercept. This accounts for between-day autocorrelation 
that might vary year to year and is unique to each radar station but is broadly 
tied to seasonality. 4) A random intercept for each radar station, accounting for 
geographic disparities unique to each station. 

Our starting models only had the first random effect without random inter- 
cepts, and we sequentially tested and retained random slopes and random inter- 
cepts based on which caused the greatest reduction in AIC. Ultimately, models for 
all response variables and seasons converged to the same random effect struc- 
ture, so we included all four random effects with random intercepts for all models. 

We then fine-tuned spline complexity by iteratively testing each fixed effect 
a) without splines and b) with three to five degrees of freedom. We chose the 
degrees of freedom for each spline based on whether a model structure resulted 
in a lowered AIC. That is, we used AIC to specify whether an included variable 
should have a spline and the degrees of freedom of this spline. Default model 
convergence criteria were changed to a gradient norm of 0.1. Spline-selection 
models often ran into boundary issues with complex random effect structures, so 
for the purposes of spline selection, we simplified the random effect structure to 
include only a nightly term unique to each radar station with a random intercept 
and a seasonal term without a random intercept. 

Once we finalized our "base" models for each season and response variable, 
we added all previously mentioned random effects and ABmax as a fixed main 
effect and with an interaction with cloud cover. In altitude models for both 
seasons and crosswind models for the fall, our models encountered singular 
boundaries in their random effects that were resolved by removing the random 
intercepts for the problematic random effects. We used log likelihood tests 
(function “anova") to test the explanatory value of models containing ABmax 
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compared to the base models. We then evaluated Cl across values of ABmax to 
assess model support for predicted trends. We used the function "sim" (package 
“arm") (37) to generate 1,000 simulated fixed effect coefficients per model. We 
used the 0.025 and 0.975 quantiles to plot the 95% CI. When model CI were 
narrow enough around a trend such that the intervals did not overlap between 
the median ABmax (~17 nT) and the extremes, we considered the model to 
show strong support for the trend. 


Gradient Boosted Decision Tree Models. We residualized our response 
variables using MLT models (function “xgb.train,” package "xgboost") (38) built 
with our weather, geographic, and temporal variables. To avoid overfitting, we 
trained a model on a train dataset (70% of data, split randomly by unique days 
in each year) until new runs failed to further the model’s explanatory power 
ona validation dataset (the remaining 30% of data) for 100 consecutive runs. 
We then fit a model on the complete dataset but restricted it to the number of 
runs identified with the validation dataset. Prior to running our final models, 
we fine-tuned the xgboost parameters by running models using a grid search 
using the training and validation sets. We kept the learning rate constant at 
0.01, but varied gamma (values at 0, 1, 3, and and tree depth (values at 6, 8, 
10, 12,14, and 16, although we restricted our final tree depths to a maximum 
of 14 to ease run times). We stopped each model after 10 rounds in which 
the model failed to improve predictions on the validation dataset. We based 
our final parameter choice for each model on the parameter combination that 
produced the lowest model R’, calculated as the covariance of the predicted 
results and the original validation dataset response variable divided by the 
variance of the original validation dataset response variable. 

To examine associations between the residualized response variable and 
ABmax, we used least square regressions with robust covariance matrix esti- 
mated confidence bands (constructed with the R package sandwich") (39, 40) 
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with the interaction between unique day and year as a grouping variable. To test 
for the effect of interactions, we excluded the interacting variable from the xgboost 
residualization and included it with a spline and an interaction with ABmax in the 
regression with robust confidence bands. Model support was assessed similarly 
as with our NLME models. 


Data, Materials, and Software Availability. Data and code for analyses are 
available at ref. 71. Code for the interpolation of ABmax is available at ref. 72. 
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